Serial No. 09/805,606 



IN THE SPECIFICATION : 

Please amend paragraph number [0024] as follows: 

[0024] Ftgr -FIG. 1 is a pictorial diagram of a computer system which comprises a 
presently preferred embodiment of the apparatus according to the invention, and upon which the 
presently preferred version of the inventive method may be carried out. 

Please amend paragraph number [0025] as follows: 

[0025] Fijgr -FIG. 2 is an illustrative example of a particle pack construction according to 
the preferred implementation, wherein spherical particles are placed along the central string. 

Please amend paragraph number [0026] as follows: 

[0026] Figr- FIG. 3 shows another illustrative example of a particle pack construction 
according to the preferred implementation, wherein spherical particles are placed along the 
central string, but wherein the particles extend more several particle diameters out from the 
central string. 

Please amend paragraph number [0027] as follows: 

[0027] Figr -FIG. 4 shows a cylindrical space with cylindrical subspaces therein 
according to the preferred embodiment and the preferred method of the invention. 

Please amend paragraph number [0028] as follows: 

[0028] Fifr -FIG. 5 shows an illustrative particle pack constructed using the preferred 
implementation. 

Please amend paragraph number [0031] as follows: 

[0031] A particl e "particle," as the term is used h e r e in herein, is used according to its 
common but broad meaning in the field. It includes any type of particle that might be found in a 
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particle rigid body of arbitrary shape containing materials. In the preferred implementation, 
particles comprise rigid or nearly rigid objects of arbitrary size and shape. 

Please amend paragraph number [0034] as follows: 

[0034] The -IUustrated in FIG. L the apparatus according to the invention preferably 
comprises a computer system 100, such as a commercially available personal computer, small 
business computer, engineering work station, mainframe, microprocessor, or virtually any other 
machine with similar processing power, appropriately configured and programmed with 
computer software or "code" as outlined and described herein for performing the tasks identified 
herein. The computer system 100 includes a processor 102, storage in the form of RAM 104, a 
hard drive 106, a compact disk (CD) drive 108, and a drive for diskettes 110, all contained within 
a cabinet 112. System 100 also includes a monitor 1 16, a keyboard 118, and a pointing 
device 120 such as a mouse, track ball or the like. System 100 also may be part of a network, in 
which case it may include a network connection 122 and network bus 124. The method 
according to the present invention, in its presently preferred implementations, may be carried out 
on computer system 100, using the same hardware and software, as described herein with respect 
to the preferred embodiment of the apparatus according to the invention. 

Please amend paragraph number [0035] as follows: 

[0035] "Machin e " "Machine," as the term is used h e rein herein, refers broadly to a 
computer, processor, microprocessor, a network of such devices, or other apparatus capable of 
performing processing as described herein. Such machines typically and preferably will 
comprise a computer, such as computer system 100. 

Please amend paragraph number [0037] as follows: 

[0037] The architecture of the preferred implementation is outlined below. In its 
presently preferred form, it is written in double precision C++ or Fortran 90. This is not, 
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however, limiting. The code may, for example, be written in FORTRAN, or other suitable 
programming languages. An example of such a computer system is shown in Fig. FIG. 1 . The 
code in its preferred version builds a single particle pack and outputs user-specified pack 
microgeometry and statistics. It is modular and can be disassembled into its component 
subroutines for use elsewhere. 

Please amend paragraph number [0090] as follows: 

[0090] The space and subspaces are constructions that define the regions or volumes 
into which the respective particles are to be placed. The space may comprise any contained 
geometric region, but preferably comprises a cylindrical construction. The subspaces also may 
comprise any contained geometric region, but preferably have a shape and configuration that 
corresponds to the overall space. In the preferred versions of the apparatus and method, the 
space comprises a cylindrical volume, and the subspaces comprise a plurality of concentric 
cylinders, disposed about a longitudinal or z axis, wherein sequential ones of the cylinders have a 
radius that is larger than the prior cylinders. An illustrative example of such a concentric 
cylinder configuration is shown in Fig. FIG. 2. The number of cylinders in the preferred 
implementation is equal to the number of particle categories, e.g., submodes, that will be present 
in the particle pack. The radius of each cylinder is selected to correspond to the radii of the 
particles. The radius of the first cylinder corresponds to or is a function of the particle radius of 
the smallest mode, the radius of the second cylinder is equal to the particle radius of the second 
smallest mode, and so on. The cylinder radius for a particular particle category preferably 
corresponds to or is a function of the radius or a representative radius of the corresponding 
particle for that mode or submode. The cylinder radius may be any positive number multiple of 
the particle radius. Preferably, but optionally, the cylinder radius is about 3 to 10 times the 
particle radius, and more preferably, about 5 to 10 times the particle radius. 
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Please amend paragraph number [00177] as follows: 

[00177] In the presently preferred implementation, a catch net is incorporated into the 
space, below which the south pole of the falling particles cannot descend. When a particle's 
south pole hits the catch net belonging to its cylinder, it is placed in the pack at that point, even if 
it is not touching any other particle. Preferably there is a separate catch net for each particle 
submode, and thus for each subspace or cylinder. Each catch net preferably covers the entire 
cross s e ction cross-section of that submode' s cylinder. If the current sphere's south pole hits 
that catch net, it is placed there. It does not see or respond to any other submode' s catch net. 
The altitude of the catch net preferably is determined dynamically, that is, it can depend on which 
placed particles are being contacted. 

Please amend paragraph number [00181] as follows: 

[00181] In the presently preferred implementation, the catch net extends across the cross 
se ctional cross-sectional area of the space, across the cross s e ctional cross-sectional areas of 
each of the subspaces. 

Please amend paragraph number [00182] as follows: 

[00182] In this preferred implementation, the catch net comprises N subnets, one of the 
subnets corresponding to each of the subspaces. Each of the concentric cylinders 1, 2, 3, ... N 
has a circula r cross s e ction cross-section p erpendicular to the central string. In the preferred 
implementation, the catch net comprises a plurality N of subnets, or individual catch nets, one for 
each of the cylinders. Each of the subnets preferably extends over the cross s e ctional 
cross-sectional area of the corresponding subspace or cylinder. The position or level of the catch 
net or subnet for each of these subspaces can and typically will differ from one another, although 
in some instances they may be at the same or nearly the same level. 
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Please amend paragraph number [00213] as follows: 

[00213] In the presently preferred implementation, the water level altitude for the z-th 
cylinder (i> 1) is defined as the average exposed north pole altitude of all particles with radius 
less than az. The user may raise or lower this level by adding a user-defined offset, but the 
default for this offset in this implementation is zero. This water level exists between the z-th 
cylinder wall of radius Wi and the next smaller cylinder wall of radius Wi.j. Thus Thus, the 
spaces between different cylinder walls have different water levels. The cylinder corresponding 
to the smallest particle's submode has no water level, only a catch net. 

Please amend paragraph number [00214] as follows: 

[00214] If a larger sphere is descending onto the pack near the center cylinder, it may hit 
the center cylinder's pack of small particles and tend to roil off its edge and fall to the floor of the 
pack or fall until it hits another large sphere. Thus Thus, the pack would segregate with larger 
particles outside smaller particles. With the water level in place, however, if the larger sphere 
rolls off the pack in the central cylinder, it simply "floats" on the water level which is at the same 
level as the "surface" of the central cylinder pack. Thus, any time a descending sphere hits a 
water level, it is placed where it is, just as was done with the catch net. Ordinarily, water levels 
will be higher than catch nets, although this is not always true. 

Please amend paragraph number [00217] as follows: 

[00217] A pack sphere, e.g., the y-th sphere, is a candidate for being the first object 
encountered by the descending current sphere if their two radii overlap, z.e., if 

(X e - Xj) 2 + (Y c - Yjf < (a c + aj) 2 ^ (35) 

If there were no other obstacles in the way, the current sphere would strike Sphere j at an altitude 
of 
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ZcO) = Zj + [{a c + aj) 2 - (Xc -Xj) 2 - (Y c - Yjff 2 . (36) 

Please amend paragraph number [00251] as follows: 
[00251] The solutions are 



-B± JB 2 -AC 

s = _ (68) 

\fB 2 -AC is negative, there is no solution. If it is zero, there is only one real solution which 
means the roll corridor just grazes the /-th cylinder, in which case it is considered a miss. If it is 
positive, there are two real solutions. The preferred implementation solves for both of them, and 
checks to determine whether either one lies in the open interval [0, 1]. If one or both lie in this 
interval, the /-th mode spheres are candidates for touching the descending current sphere. 

Please amend paragraph number [00267] as follows: 

[00267] The preferred implementation also can find which water level the current sphere 
might hit. Water level altitudes between adjacent cylinder walls are generally different as the 
pack is constructed. The only condition under which the current sphere might contact the £-th 
water level Z w (k) is if and only if the distance P c of the current sphere from the global Z axis lies 
between cylinder wall radii Wk-i and Wk at the point of contact. P c is given by 

P C = (X C 2 + Y c 2 ) m (83) 

where X c and Y c can be obtained from Eq. 80 as 

X c = Xi + (a c + a { -)(-cos0y/ cos^- sin# c ' cos#> c ' + sin^ sin# c ' sin^' 
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+ sin#; 



)i cosfpji cosOc) 



(84a) 



and 



Y c = Yi + (a c + a/)(-cos# y 7 sin^y/ sin# c ' cos^ c ' - cos^ sin# c ' sin#> c ' 
+ sinOji sirups cosd c ') . 



(84b) 



Please amend paragraph number [00287] as follows: 

[00287] The square root has a positive sign because the difference between q> c and cp u can 
be no greater in magnitude than n/2. Again, the arccosine has two solutions: the principal 
value Cos" 1 and its negative -Cos" 1 . The current sphere would lose contact with the upper sphere 
if it is at azimuthal position cp u ± Cos" 1 . Its true position depends on where it was initially. If its 
initial position is (p C j > (p U9 then the principal value Cos" 1 is used. If its initial positionjs 
<Pci < <Pu> then the negative of the principal value is used -Cos" 1 . In the unlikely case that the 
initial position equals q> u , then the current sphere is balanced on the upper and lower sphere and a 
direction is randomly chosen by the random number generator. The preferred implementation 
using these relationships is to determine where contact is lost. 

Please amend paragraph number [00298] as follows: 
[00298] Squaring both sides gives 



(A 2 - B 2 ) 2 sinV - 2Bj(A 2 - B 2 ) sinV + (A } 2 - 2AoA 2 - 2A 2 + B } 2 + 2AoB 2 + 
2A 2 B 2 ) sin V + 1(Ao + A 2 )Bj sin^ c ' + (A 0 'Ai+A 2 )(A 0 + Aj+A 2 ) = 0^ 



(110) 



Please amend paragraph number [00314] as follows: 
[00314] This equation is of the form 

A 0 +Aj oos<pc + A 2 cos (pc + Bj sirups + B 2 sin <p c ' = 0 



(132) 
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where in this case, the coefficients are given by 

Aq^Xi + Y/ 2 + 2(a c + a t ) cos# c '(^G sin0 y/ cos^ yi - + Y ( sin^y/ sin^) + 



(a c + a)) 2 cos 2 9 c ' sin 2 0 y/ (133a) 

-4/ = -2(a c + A,-) cos^/ sin# c ' cos^-,- + 7/ sin^-,- + 

(a c + a,-) cos# c ' sin^,/] (133b) 

^ = (a c + a/) 2 cos 2 0 y/ sin 2 0 c ' (133c) 

5y - 2(a c + a,) sin# c ' {sintf,, [X t + (a c + a/) cos# c ' sin# y/ cos^-J 

- cos# c 'jY/ + (a c + a,-) cos# c ' sin# y/ sin^-,-]} (133d) 

B 2 = (a c + fl/) 2 sin 2 0 c f ^ (133e) 



Please amend paragraph number [00326] as follows: 

[00326] In this expression, W c is the radius of the cylinder wall, P t = ^JX? + Yf is the 

radial position of Sphere / from the global Z axis, and a c and a, are the radii of the current sphere 
and Sphere /, respectively. The plus sign is used in Eq. 138 if the current sphere is rolling toward 
increasing 0, the minus sign is used if the current sphere is rolling toward decreasing 0. If the 
absolute value of the quantity in the parentheses of Eq. 139 is greater than unity, then there is no 
solution and the current sphere simply rolls around the perimeter of the cylinder until it reaches 
the lowest point it can. This lowest point is opposite the azimuthal position of Sphere L H e nc e 
Hence, if 
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w c 2 + p? -o 



(a c + a,) 2 



>1 



(140) 



2KP, 



then 



A0 c -±n . 



(141) 



Please amend paragraph number [00330] as follows: 

[00330] The current sphere will impact a placed pack sphere, say the y'-th sphere, if there 
is a solution to 



anywhere in the range (0 c o, 0 C f)- The Cartesian components of the vector /? c (O c ) = X c , Y c , Z c are 
given by 



Please amend paragraph number [00345] as follows: 

[00345] When the pack is still small enough that an insufficient number of particles have 
been placed to satisfactorily define a surface, the catch net may be ratcheted up a little after the 
placement of each particle so that the particles on the bottom of the cylinders do not all have the 
same altitude. If the switch ROUGHNET is true, then the current sphere's catch net moves 
upward by a small random amount 



fl e (O e ) - R\ = a c + aj 



(142) 



X c (0 c ) = W c cos0 c 



(143) 



Y c (® c ) = W c sm0 c ^_ 



(144) 
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ra min IN cs (157) 

where r is a uniformly generated random number in [0, 1] and N cs is the total number of sphere 
cross s e ctions cross-sections required to cover all the cylinders: 



tf fl = S Wk f*' 1 (158) 



where n is the number of modes, Wk is the cylinder wall of the &-th mode, and Wo = 0. 
If ROUGHNET is false, then the catch net stays where it is until enough particles have been 
defined to make a surface and until the surface is more than 2 sur f ace - yak above the starting place. 

Please amend paragraph number [00362] as follows: 
[00362] The sphere is completely outside the cylinder 

V in = 0 (162) 

only if one of the following conditions are met: 

If the sphere's north pole is lower than the cylinder bottom . . . 

z s + a s <zi> (163a) 

Qf4£ or if t he sphere's south pole is above the cylinder top . . . 

z s -a s >z t (163b) 

Qptf or if the sphere is more than its radius away from the cylinder side . . . 
p s -a s >a c (163c) 
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Qf4f or if the sphere center is above the cylinder top and is more than its radius 

away from the upper cylinder edge. . . 
(p s - «c) 2 + fe - ztf > a* 2 and z s > z t (163d) 

Or if or if the sphere center is below the cylinder bottom and is more than its 

radius away from the lower cylinder edge . . . 
(p s - «c) 2 + (zt - z s f > a s 2 and z s <z b . (163e) 



Please amend paragraph number [00363] as follows: 
[00363] The sphere is completely inside the cylinder 

Yin = (4/3)W (164) 

only if all three of the following conditions are simultaneously met: 

If the sphere's north pole is no higher than the cylinder top . . . 

z 5 + a 5 <z t (165a) 

And and if the sphere's south pole is no lower than the cylinder bottom . . . 
z s -a s >z b (165b) 



And and if the sphere's equator nowhere exceeds the cylinder radius . . . 

p s + a s <a c . (165c) 

Please amend paragraph number [00364] as follows: 
[00364] The cylinder is completely inside the sphere 

Vin = na c 2 (z t - z b ) (166) 
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if both the following two conditions are met: 

tf-fee -if the upper edge of the cylinder is inside the sphere . . . 

(ps + a c f + (z, - z s f < a 2 (167a) 

And and if the lower edge of the cylinder is inside the sphere . . . 

(ps + a c f + (z s - z b f < a 2 . (1 67b) 

Please amend paragraph number [00365] as follows: 

[00365] There are several subsets of Case 4. To efficiently organize them, we define 
two circles: the sphere circle and the cylinder circle. The sphere circle is the perimeter of the 
sphere's cross s e ction cross-section at a given z altitude. Its radius is 

r s = [a s 2 - (z - z s f} m , z s - a s <z<z s + a s (168) 

and its radial position p s is given by Eq. 161. Likewis e Likewise, the cylinder circle is the 
perimeter of the cylinder's cross s e ction cross-section which has constant radius a c at anyz 
altitude between and z t and is centered on the z axis. 

Please amend paragraph number [00388] as follows: 
[00388] Lastly, 

E = A/3 [{A 2 ^B 2 )E{X h q)-2B 2 V(k h q)] - f 7 /3 [(A 2 - tfygi 2 • B 2 )] m - 

A/3 [(A 2 + B 2 ) E(l 2 ,q) - IB 2 F(^)] - fc/3 [(A 2 - ^{Ci - B 2 )f\ £ 2 <A^ (202) 

Again, the entire second line of this expression is zero if Q = A. 
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Please amend paragraph number [00394] as follows: 
[00394] The sphere is completely outside the cylinder 

V in = 0 (212) 

if one of the following conditions are met: 

If the sphere's north pole is lower than the cylinder bottom . . . 

z s + a s <z b (213a) 

Or if or if t he sphere's south pole is above the cylinder top . . . 

z s -a s >z t (213b) 

Qi4£or if the sphere is more than its radius away from the cylinder side . . . 

p s -a s >a c (213c) 

Or if or if the sphere center is above the cylinder top and is more than its radius 
away from the upper cylinder edge . . . 

(p s - «c) 2 + (z s - z t ) 2 > a 2 and z s >z t and p s >a c (213d) 

Or if or if the sphere center is below the cylinder bottom and is more than its 
radius away from the lower cylinder edge . . . 

(p s -ac) 2 + (z b -z s ) 2 >a s 2 and z s <z b and p s >a c . (213e) 
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Please amend paragraph number [00395] as follows: 
[00395] The sphere is completely inside the cylinder 



^ = (4/3)W (214) 



only if all three of the following conditions are simultaneously met: 

If4ke -if the sphere's north pole is no higher than the cylinder top . . . 

z s + a s <z t (215a) 

And if and if the sphere's south pole is no lower than the cylinder bottom . . . 

z s -a s >z b (215b) 

Attd-tf and if the sphere's equator nowhere exceeds the cylinder radius . . . 

p s + a s <a c . (215c) 

Please amend paragraph number [00396] as follows: 
[00396] The cylinder is completely inside the sphere 



V in = %a c \z t -z b ) (216) 



if both the following two conditions are met: 

If-fehe- if the upper edge of the cylinder is inside the sphere 
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(p s + a c ) 2 + (z t -z s ) 2 <a s 2 



(217a) 



Aftd4f-and if the lower edge of the cylinder is inside the sphere . . . 



(p s + a c ) 2 + (z s -z b ) 2 <a s 2 . 



(217b) 



Please amend paragraph number [00397] as follows: 

[00397] Again, there are several subsets of Case 4. To efficiently organize them, we 
define two circles: the sphere circle and the cylinder circle. The sphere circle is the perimeter of 



and its radial position p s is given by Eq. 211. Lik e wis e Likewise, the cylinder circle is the 
perimeter of the cylinder's cross s e ction cross-section w hich has constant radius a c at anyz 
altitude between z& and z, and is centered on the z axis. 

Please amend paragraph number [00402] as follows: 

[00402] Armed with these four solutions, we will now consider each of the four subcases 
of Case 4, but we will consider separately whether the sphere circle lies in the southern 
hemisphere of the sphere or whether it is the equatorial plane or northern hemisphere. We 
consider the four subcases separately for the two hemispheres making eight situations in all. We 
consider the hemispheres separately because if the sphere circle is in the southern hemisphere, it 
is increasing in size as z increases. -Henee - Hence, the preferred implementation can determine 
the altitude at which subcases change. If the sphere circle is its equatorial plane or is its northern 



the sphere 's- 



cross-section at a given z altitude. Its radius is 



r s = [a* 2 ~(z- z s f] m , z s -a s <z<z s + a s 



(218) 
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hemisphere, it is decreasing in size as z decreases, and the altitude at which subcases will change 
again can be determined. 
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